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Abstract 

We present a new method of obtaining time-dependent matrix elements of the 
■ electromagnetic pulse produced by a highly-relativistic projectile. These matrix 

elements are used in a coupled-channel calculation to predict the cross-sections for 
i-H \ population of 1- and 2-phonon states of the giant dipole resonance. Comparisons 

are made with the predictions of the long- wavelength and Born approximations. 

a \ 

: 1 Introduction 

The subject of relativistic Coulomb excitation has received extensive study in the two 
decades since its theoretical foundations were established in the classic work of Alder and 
Winther pQ. In the semiclassical approach to relativistic nuclear Coulomb excitation, the 
relative motion of the projectile and target is treated classically. Indeed, it is usually 
assumed to be straight-line motion at constant speed v. The evolution of the internal 
degrees of freedom of each nucleus, under the influence of the classical electromagnetic 
field produced by the other nucleus, is then followed using quantum mechanics. 

There exists a wide literature devoted to the use of relativistic Coulomb excitation 
for the study of several aspects of nuclear structure that cannot be explored through the 
nuclear interaction [2-8]. In recent years, this subject has been also shown to be relevant 
to the practical question of the stability of beams of relativistic heavy ions, since some of 
the processes which lead to loss of beam ions are initiated by Coulomb excitation [2,9-11]. 

The situations that are easiest to interpret are those that can be described in terms 
of the Coulomb excitation of a single level. If this level is unbound, the Coulomb ex- 
citation will be followed by particle emission. Methods have been developed to extract 
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from the particle emission cross-section the part that is associated with the Coulomb ex- 
citation of the emitting state, as opposed to excitation via nuclear forces [121 ESI- The 
Fermi- Williams- Weizsacker (FWW) method of virtual quanta ^1] is frequently used to 
analyze processes in which Coulomb excitation is followed by particle emission (see, e.g., 
refs. [5,12,15,16]). Here one simulates the electromagnetic pulse by an equivalent flux of 
virtual photons, and uses experimentally determined photo-nuclear cross-sections to de- 
scribe the effect of these virtual photons on the target. However, most of the FWW 
applications to relativistic Coulomb excitation at energies of spectoscopic interest have 
not used experimentally determined photo-nuclear cross- sections, but have used calcu- 
lated electromagnetic matrix elements. Thus, in effect, these calculations really amount 
to the use of the Born approximation as developed in pQ, although they are expressed in 
the language of the method of virtual photons. 

Analyses of experimental data based on the Born approximation imply that the final 
state is weakly coupled to all states other than to the initial state. The Born approximation 
is just the first term of the Born (or Fredholm) series, which constitutes a formal solution 
of the time-dependent Schrodinger equation. In dealing with multi-channel processes, 
some workers J7| have attempted to include a few higher terms of this series. Other 
approaches, based on coupled-channel methods, [SJE], have attempted a numerical solution 
of the Schrodinger equation which, in principle, sums all terms of the Born series. The 
harmonic vibrator method allows for complete coupling of all the states of a vibrational 
band, if some simplifying assumptions are made concerning the ratios of interaction matrix 
elements. 

All these calculations require matrix elements of the electromagnetic interaction be- 
tween the nuclei, defined with respect to the eigenstates of internal motion. If the relative 
velocity of the nuclei is relativistic, the interaction is strongly retarded, which introduces 
difficulties into the evaluation of its matrix elements. The long- wavelength approximation 
(LWLA) has proven to be a very useful device for overcoming these difficulties jSJ E|. 

We see then that two approximations have been extensively used in this field: the Born 
approximation to solve the Schrodinger equation and the long-wavelength approximation 
to evaluate matrix elements. These approximations have been found to be effective below 
bombarding energies of 2 GeV per nucleon, which has been adequate for most nuclear 
structure investigations done so far. At higher energies, which could be of interest in 
future studies, and for energies presently considered in colliders, these approximations 
may not be adequate. 

Much of the recent effort in the theoretical study of relativistic Coulomb excitation 
is devoted to the nuclear structure aspects of the problem. For example, in the study 
of giant resonances, one must take account of the spreading of the resonance amongst 
the background states (see, e.g., refs. [SI IHj ITHl IT§1 HU[ l2*Tj). and the anharmonicity of the 
oscillation [H [HI [221 12B1 121] • Of course, these studies also must incorporate the interaction 
of the nuclear motion with the electromagnetic pulse, and this is generally done by means of 
the Born approximataion and/or the the LWLA. For this reason, we believe it is important 
to understand the limitations of these approximations, even in Coulomb excitation studies 
whose primary emphasis is on nuclear structure. 

The purpose of this article is to assess the ranges of validity of the LWLA and the Born 
approximation. We accomplish this by comparison with a general numerical solution 
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of the time-dependent Schrodinger equation, expressed as a set of coupled differential 
equations. The energy range considered here extends from energies currently used for 
nuclear structure studies (see, e.g. refs. [5,6]), up to about 10 GeV per nucleon. The 
methods employed are quite general, and could be used at still higher energy. 

In Section 2, a general survey of methods and approximations is given. In Section 3, we 
present our method for numerical evaluation of the Fourier transform of interaction matrix 
elements. In Sections 4 and 5, we illustrate our procedure by investigating the relativistic 
Coulomb excitation of giant dipole resonance (GDR) phonons. We present numerical 
results for the particular case of 208 Pb projectiles exciting GDR states of 40 Ca. These 
results are compared with those produced by approximate approaches. We summarize in 
Section 6 the main results we obtained. 



2 Survey of methods used in RCE 

In the semiclassical theory of relativistic Coulomb excitation, the electromagnetic field 
between the projectile and target is treated classically, but their internal degrees of free- 
dom are treated according to the principles of quantum mechanics. The time dependent 
Schroedinger equation is written 

ih ?M = [ Ho + V(t)]ij(t) (2.1) 

with Hq referring only to the target 1 degrees of freedom, and V(t) the interaction between 
the target and the electromagnetic field of the projectile: 

V{t) = J [ V f(v\t)p(r>) - - c Af(r',t)-i(r')]d 3 r' . (2.2) 

The scalar and the vector potentials associated with the projectile electromagnetic field are 
ifiQ 1 ^, t), Ag; t (r',t), and p(r'), j(r')are the target transition charge and current densities, 
respectively. 

The projectile is assumed to have a spherically symmetric charge distribution, with 
total charge Zpe. Because of its large momentum, it can be assumed to follow a straight- 
line trajectory at constant speed v. Thus its center is located at time t by 

r = by + vtz . (2.3) 

The constant b is the impact parameter. Excitation probabilities are calculated as func- 
tions of b. The scalar and vector potentials of the projectile field, at the target point r', 
are given by jHj 



yV 2 + (y'- b) 2 + 7 2 (^' - vt)- 



Af(r',t) = V c ct (r',t)z (2.4) 



■"■For simplicity, we restrict our discussion to the situation in which only the target is excited, but the 
argument is easily generalized to allow for projectile excitation as well. 
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Because of the factor 7 2 in the denominator of Eq.(2.4), these potentials are not spherically 
symmetric about the projectile center. 

The Schroedinger equation (2.1) is conveniently expressed in terms of an expansion in 
eigenstates <p a of the nuclear Hamiltonian H , 

m = J2 Oa(*)e* e °Va , (2-5) 

a 

with 

H (f)a = £ a <Pa , < 4>a\4>f3 > = <W J (2.6) 

Vap{t) = < a |/ c et (t)p - -A r c et (t)-j|<^> . (2.7) 

c 

If (2.5) is substituted into (2.1), the result is a set of coupled ordinary differential equations 
for the amplitudes a a {t)\ 

. h da$l = ^ e , (ep _ ea)t y pa[t)aa{t) (2 8) 

These equations must be solved subject to initial conditions 

a a (-oo) = 5 afi . 

The probability that the target will be in state 0^ after the collision is given by 
|a j g(+oo) | 2 , and the cross-section for the population of 4>p is given by 

poo 

op = 27r|a /3 (+oo)| 2 M& . (2.9) 

" ^min 

Here 6 min is usually taken to be somewhat larger than the sum of the target and projectile 
radii [3. 

The first step in the solution of Eq.(2.8) is the computation of the interaction matrix 
elements Vp a (t). This is a formidable task, since the retardation of the interaction intro- 
duces a directional asymmetry, which means that the interaction is not invariant under 
rotations of the internal coordinates alone. Fortunately, Alder and Winther have found a 
convenient multipole expansion of the Fourier transform of the Vp a (t), 

Vp a (u) = / -e^% a {t) . (2.10) 



The Vp a {u)) defined in Eq.(2.10) can be used directly in a coupled integral equation for- 
mulation of the Schroedinger equation, as was done in ref.[22]- However this approach 
is difficult to implement when many states cj) a have to be included in the calculation. 
Another alternative is to proceed by inverting the Fourier transform (2.10) 

V 0a (t) = / due-^% a (u) (2.11) 

to convert the Va a (o;) obtained from the Alder- Winther expansion into Vp a (t) which can 
be used in the time- dependent formulation (2.8) of the Schroedinger equation. 
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There is no known closed expression for the Fourier transform (2.11). In ref. [H] Lanza 
et al replaced the exact Alder- Winther expression for Vp a (uS) by its long- wavelength ap- 
proximation (LWLA). Here one assumes that the important Fourier components of the 
electromagnetic pulse correspond to wavelengths that are large compared to target di- 
mensions. In this case, the small-argument limit of the spherical Bessel function, 



can be used. With this approximation, Lanza et al [S| were able to obtain explicit ap- 
proximate expressions, in terms of hypergeometric functions, for the Vp a (t), which they 
used in their analyses of multiphonon Coulomb excitation. They tested the validity of the 
LWLA by comparing a few Vp a (t) calculated using the LWLA with numerical evaluations 
of the Fourier transform (2.11) of the exact Vp a (u). They concluded that the LWLA was 
adequate for their analysis of the 208 Pb+ 208 Pb collision at Ei a b=641 MeV. 

Bertulani et al [3] approached the problem of finding approximate Vp a (t) by expressing 
the potential (2.4) in a Taylor expansion around x' = y' = z' = 0. This formally exact 
expansion was truncated in a manner that caused the remaining terms to be precisely equal 
to those given by the LWLA expressions in ref . jB] . However, the advantage of the truncated 
Taylor expansion is that it yields simple expressions for the Vp a (t) which, although they 
are equal to the hypergeometric functions used in ref.jHj, are much more transparent. The 
general proof of the equivalence of these methods can be found in ref. [2^] • 

The approach to be followed in this manuscript is the use of a quick and accurate 
method for the numerical evaluation of the Fourier tranform (2.11), using exact Vp a (u). 
The Vp a (t) calculated in this way will be used in the numerical solution of the coupled 
time-dependent equations (2.8). We will thus be able to obtain essentially exact Coulomb- 
excitation cross-sections at any bombarding energy. These exact cross-sections can be used 
to explore the limits of validity of the LWLA and other approximations. 



3 Numerical Evaluation of the Fourier transform 

The difficulty of evaluating (2.11) as a numerical integral is the rapid oscillation of the 
integrand for high values of t. Our approach to this problem is a generalization of the 
idea behind the use of Simpson's rule in the evaluation of f* f f(x)dx. We first divide the 
cu-integration range of (2.11) into an even number of steps of length h. We assume that 
these steps are small enough so that, over every adjacent pair of steps, Vp a (u) can be 
approximated by a quadratic function of uj. Thus for the interval U\ — h < uj < U\ + h, 
we make the approximation 

Vp a {uj) ~ (— — ^ 2 u + {—r-^) v +w 
n n 

We can make this approximation exact at u = u\, uu\ ± h by defining u, v, and w to be 

Vp a {ui-h) + Vp a {ui + h) 

U = VpaiUi) 
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VpJfJi + h) - Vpgjui - h) 
V 2 

W = Vpa{u)i) ■ 

Note that since the quadratic approximation is applied to Vp a {uj), its validity is indepen- 
dent of t. Then the uji — h < uj < ui + h part of the integration (2.11) can be approximated 
by 



due-^% a (u) ~ r +h due-^ [ it^ 1 ) 2 u + (^ T ^ i ) v + w } (3.1) 



/ due- luJ % a {u) ~ (3.2) 



rwi+h rui+h 

due- tuJ % a (iu) ~ / 

lu)i—h Juii—h 

Now an exact integration of the right-hand side of (3.1) yields 

(■uii+h 

2 



1 'fs{ htcos(ht)(2ucos(uJit) + vht sm.(uit) ) 
+ sm(ht) ( (—2m + h 2 t 2 (u + u>)) cos(c<Jit) — i>M sin(c<Jit) ) | 
-| ht cos(ht) (2u sm(uit) — v ht cos(c<j 1 t) ) 



+ sm(ht) ( (— 2-u + h 2 t 2 (u + tu)) sin(c^it) + vht cos(cj]t) ) | 

This approximation is used for every adjacent pair of steps in the uj integration. 

Equation (3.2) cannot be used at t=0, since it takes the form 0/0 there. By expanding 
the quantities between { } about t—0, we can show that 



ru>l+h 

lim / due~ l,JJ % a {u) 



2 t 2 
~ -h(u + 3w) - — (h 2 (3u + 5w) + lOhvuJx + 5(it + 3u)w?) 

— ^/it( /iu + (u + 3iu)wi) (3.3) 

The only circumstance where the quadratic representation of Vp a (u) is inadequate 
occurs when Vp a (u) is singular at u=0. In these cases Vg a (u;) has the form 

^» = f(u) Jf (l^) , (3.4) 



7f 



with /(0) 7^ 0. Since K (— ) diverges logarithmically as c<j — > 0, numerical integration of 



(3.4) requires special precautions. In this case, we work with the identity 

Vpaiw) = (f{u) - /(0)) Kofe^) + f(0) K (^) . (3.5) 
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The first term in (3.5) is regular at u = 0, and its Fourier transform can be evaluated 
without difficulty using Eqs.(3.2) and (3.3). The Fourier transform of the second term in 
(3.5) can be evaluated exactly 



This term dominates the interaction matrix element Vp a (t) at large values of \t\. 

4 Brink's model for the Giant Dipole Resonance 

4.1 The GDR phonon states 

The giant dipole resonance can be regarded as a collective oscillation of the protons in a 
nucleus relative to the neutrons [21] ■ The one-phonon GDR state also has a simple inter- 
pretation in terms of an isovector linear combination of one-particle, one-hole excitations 
of the ground state [2Ej- It was shown by Brink j2Hj that states of the GDR have a simple 
interpretaion within the harmonic oscillator shell model, in terms of harmonic oscillations 
of the vector R pn connecting the centers of unexcited proton and neutron spheres. We 
use this representation of the GDR to provide the transition charge and current densities 
needed for the calculation of the Vp a {uj) (see App. A). 

An individual GDR state is described by the three quantum numbers N, L, M, in which 
N(= 0, 1, 2 • • •) is the principal quantum number. The number of GDR phonons associated 
with this state is 2N + L. Thus the j3, a in Vp a {uj) stand for Np, Lp, Mp and N a , L a , M a , 
respectively. These states are denoted ^m L - 

4.2 Structure of the matrix elements 

The algebraic details of the construction of the GDR states, and the derivations of the 
formulae for the associated Vp a (u), are given in the Appendices. In order to facilitate 
our discussion of the long- wavelength approximation, and of the features that lead to 
convergence of the uj integral, we summarize here the general structure of the result. 

Vp a (uj) can be expressed in terms of the Fourier transforms of the scalar and vector 
potentials as follows: 



Vp a (uj)= <*Z P \ V ^M a La > = /rfV[p, a (r0^ ct (r / ,^)--j, a (r')-A- t (r',^)] , 

(4.1) 




(3.6) 



where 




The multipole expansion of (p™ 
written in the form 




<f (r» = E *m(^!?) Ca.» iA(^r') K*(f) (4.3) 
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with 



where 



= 'h' A " X ^ <2A + 1} <A " <A + ? (2-,r(, i + ,).:.(A-,-2 n) . 

It is shown in Appendix B that the / </?p part of the matrix element can be expressed 
in the form: 



J d\< <tf{r> ,u)PPaV) = £^>»^) / dV P, a (r') Jx (K')Y^f') (4.4) 



\u\b 



~ e*^ £ 2(2* + 1) r jo(-r')ul, e (ry 2 dr> £ A^r ^(™)u*l(^V2), 

(4.5) 

where A n £NX are suitable coefficients. The functions u n fX r ) are the harmonic oscillator ra- 
dial functions associated with the individual nucleons moving in the shell-model potential, 



u„i(r) = 



2^n+2( 2 £ + 2n+l)!!z/ 3 /2 ^ (-2f n\ . 2 . K 

\ ^ ( ^ )e 2 *? d(n-.)!(2f + 2. + l)!! ^ > ' 

(4-6) 

where z/ = muj sm /h, with u; sm representing the frequency associated with the shell-model 
harmonic oscillator potential. The tilde in the corresponding Fourier transformed states 
u signifies that the size parameter is taken to be 2/(Ztv)- The matrix elements of the 
~ \ fjpa ■ A™* part of the electromagnetic interaction are easily represented in terms the 
above formulae (see Appendix C). 

The advantage of this approach is that it yields explicit functions for the matrix ele- 
ments Vp a (iv), with only the single parameter v (which is related to the nuclear radius). 
Of course, the representation of GDR states in terms of eigenstates of the independent 
particle shell model with oscillator radial functions is highly schematic. Some important 
features of the real GDR, such as the spreading of the GDR among background states, 
are missing from this model. However, the states we use give a realistic picture of the 
collective oscillation of the nuclear protons and neutrons relative to each other, and we 
believe they provide sufficiently accurate Vp a {uj) to enable us to test the dynamics of the 
Coulomb excitation process. 

4.3 The Long- Wavelength Approximation (LWLA) 

4.3.1 Expression of the LWLA in Brink's model of the GDR 

We can get the result of using the long wavelength approximation (2.12) in (4.4) by keeping 
only the lowest powers of ^ in (4.5), apart from the to dependence of K^M P -M a \(^~)- This 
implies the following replacements in (4.5): 
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r U^r')ulAry 2 dr' — f°° < £ (r>'W 



a; /c s 



«jvl(-t-^) 



2 L - JV + 2 (2L + 2AT+ 1)!! 



a; /c 



2^' ^(Zz//2)3/2 0F[(2L + l)!!] 2 ^v^/ 

If these replacements are made, our calculated results are equal to those obtained with 
the LWLA. In Section 5 we will test the adequacy of the long-wavelength approximation, 
for varying bombarding energies and impact parameters. 



4.3.2 General conditions for the validity of the LWLA 

If the integral representation (2.11) of Vp a (t) is to converge, it is necessary that Vp a (u) 
should decrease sufficiently rapidly as — > oo. We use the term damping to refer to 
the phenomena responsible for this decrease. It is evident from (4.4) that the decrease of 
Vp a {uj) with increasing \uo\ has two causes: 

1) Kinematic damping. The factor K^-^) has a large- \u\ limit of 



7T 

1 — rr— e f v 

iv ' 

This decrease occurs because the time-width of the electromagnetic pulse at the target is of 
the order of (^), which contains frequency components up to, but not greatly exceeding, 
w max ~ It can also be said that the impulse is too adiabatic to contain frequency 
components with u > u; max . 

2) Dynamic damping. If ^R » 1, the factor j\(^-r') in (4.4) will go through many 
oscillations as r' is integrated from to R. These oscillations will lead to cancellations 
between different parts of the r' integration range, which will lead to a decrease of Vp a (cu) 
when \uj\ > c/R. This phenomena can also be called retardation damping, because it is 
associated with the finite speed, c, of propagation of electromagnetic signals. Because of 
this finite speed, different parts of the target nucleus, at one time t, are influenced by the 
projectile at different points along its orbit. This implies that the r' integration needed 
for the evaluation of Vp a {uj) effectively produces a time average over the projectile history, 
which smooths out sharper features of this history. Thus Vp a (t) will not vary as rapidly 
with t as it would if there had been no retardation. Equivalently, the presence of high- 1 | 
components in Vp a (ou) is diminished. 

It is worth remarking that the criterion for dynamic damping, \u\ > c/R, is opposite 
to the validity condition for the LWLA, |u;| < c/R. Thus if we remain within the regime 
of validity of the LWLA, we will not experience dynamic damping; in this regime, the only 
phenomenon that can lead to the damping required for the convergence of the integral 
(2.11) is kinematic damping. Hence we conclude that the LWLA can only be applied if 
| a; | < c/R for the entire | a; | -range from u = up to c<j max ~ the frequency at which 
kinematic damping becomes effective. This implies that 

7f c 

<^max ~ -T < S > 



7f 
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b 7t> 



R c 



must be true if the LWLA is to be trusted. 

5 Numerical results in 40 Ca— 208 Pb collisions 
5.1 Comparison of excitation amplitudes 

The coupled differential equations (2.8) for the ap(t) were numerically integrated using 
the fourth-order Runge-Kutta method. We have included all states with < two phonons, 
i.e. the states with (NLM) = (000), (010), (011), (020), (021), (022), (100) (M — 1,2 here 
represents the reflection-symmetric combinations of the M = ±1,±2 states.) A test 
of the numerical accuracy of the integration procedure is the extent to which the total 
normalization of the state, \ a /3(t)\ 2 — 1> is preserved. It was not difficult to satisfy this 
condition to 8 decimal places, even in situations when there was strong coupling between 
the different ap(t). 

The only free parameter in our model is the size parameter, v = muj sm /h, that char- 
acterizes the shell-model potential. For simplicity, we follow the usual prescrition 



which, in the case of 40 Ca, leads to hu sm =11.7 MeV, v = .28177 far 2 . 

We illustrate the criterion of validity of the LWLA with the example of 208 Pb projectiles 
of kinetic energy 10A GeV, on a 40 Ca target. Then 7 = 1 + 10/. 938 ~ 11.66, and the 
condition of validity of LWL is 



Figures la and lb show Vp a {uS) and Vp a (t) for b = 12 fm, with a referring to the 
ground state and (3 referring to the one-phonon state with M — 1. It is evident that the 
exact expression for Vp a (u), which includes dynamic as well as kinematic damping, has 
less high- 1 u; I content than does the LWLA expression. Correspondingly, Vp a (t) calculated 
with the LWLA is more sharply peaked than is the exact Vp a (t). Figures 2a and 2b show 
the same comparison, but for b = 50 fm, where the above argument suggests that the 
LWLA should be adequate. While it is clear that the LWLA does a better job at b = 50 
fm than at b = 12 fm, the approximation at b = 50 fm is still only fair. 

The calculation of the total excitation probability requires an integration over b. Thus 
the inadequacy of the LWLA for small b will have serious consequences only if the small-6 
part of the integration range makes an important contribution to the total 6-integration. 
We will see in Section 5 that the bombarding energy determines whether or not this occurs. 

Some typical plots of Pp{t) = \ap(t)\ 2 are shown in Figures 3 through 6. 

Figure 3 shows the occupation probability of vl^ 1 for an impact parameter of 200 fm. 
The horizontal line indicates the Born approximation prediction. It is seen that, beginning 
near t=0 fm/c, the occupation probability rises almost monotonically to the Born value, 
and then remains there. The occupation probability of ^/[j 1 shown in Figure 4 has a very 
different behavior. The coupling between and ^q 1 has a long range, because the 



sm 



40A"3 MeV, 
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A = jj, = term in the j ■ A part of the interaction can connect these states. It is seen 
in Figure 4 that the occupation probability rises to a maximum near t — fm/c, and 
then strongly decreases. The asymptotic probability in this case is 0.6 x 1CT 5 , compared 
to the Born prediction of ~ 1CT 8 . The rise and fall of the ^(j 1 occupation probability is 
associated with the fact that ^q 1 represents a longitudinal oscillation (with reference to 
the direction of the projectile motion), and the coupling to the ground state changes sign 
during the encounter, ty® 1 , on the other hand, is a transverse oscillation, and the coupling 
to the ground state has a constant sign. 

The Born prediction for the excitation probablility of M = states is always very 
small at high bombarding energy (see Section 5.3.1 below). 

The corresponding curves for a grazing collison (b = 12 fm) are shown in Figures 5 and 
6. Again the ^Z? 1 occupation probability rises almost monotonically to a constant value, 
but in this case the constant value is significantly less than the Born probability. The 
population of ^{j 1 shows a rather complicated time structure. The asymptotic occupation 
probability is again much greater than the Born probability. 

The small-amplitude oscillation of the occupation probability of ^[j 1 shown in Figure 
6 for b = 12 fm and t > 100 fm/c is a consequence of the long range of the A = \x = 
coupling term. 



5.2 Evaluation of the cross section 

Figures 7 and 8 show the asymptotic probabilites for ^{j 1 and ^J 1 , plotted as functions of 
impact paramater b. 

In the case, it is seen that the Born expression gives a very good representation 
of the exact probability when b ~ 100 fm. In this situation, the most convenient way to 
evaluate the cross-section is to break the impact-parameter integral (2.9) into two terms: 

rb' roo 

up = I 2n\af j (+oo)\ 2 bdb + / 27r|a /3 (+oo)| 2 M&, 

J6 min Jb' 

where b' ~ 100 fm. The first integral is evaluated numerically using |a/3(+oo) | 2 calculated 
by integrating the Schrodinger equation for a range of b values between b min and b'. The 
second integral can be calculated exactly, because the b dependence of the Born probability 
is given by (f^(^)) 2 , and one can take advantage of the exact formula 

j~{K„{x)Yxdx = | [ (iWO) 2 - (^(O) 2 - yiW0^(0 ] 
Thus the contribution to the cross-section from b' < b < oo will be 

° = *q 2 b' 2 1 [ (iWO) 2 - (MO) 2 - y^WO^MO ] 

where 

_ fr'^on-shcll 

7f 
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and 

2 _ Born probability at b' 

q = OMO? 

In a situation such as that shown in Figure 8, where the Born formula does not give 
an adequate approximation to the exact result even at large b, it is necessary to do a 
numerical integration of (2.9) over the entire b range. 

5.3 Comparison of cross sections 

We use the term "exact" to refer to calculations in which the interaction matrix elements 
are evaluated without approximation, and full account is taken of coupling between 0-, 1-, 
and 2-phonon states. The LWLA calculations also include full coupling between 0-, 1-, 
and 2- phonon states, but the interaction matrix elements are evaluated approximately 
(Equation 2.12, and Section 4.3). The Born approximation uses the correct on-shell inter- 
action matrix elements, but assumes that the transition to any final state is accomplished 
in a direct, single- step process. Of course, our "exact" calculation is not really exact, 
since it leaves out the effect of coupling with states of 3, 4, ... phonons. For this reason, 
we have chosen to restrict our attention to bombarding energies of < 10 A GeV, with the 
expectation that below this energy a description involving only 0-, 1-, and 2-phonons will 
be adequate. This is probably a fairly good approximation for the excitation cross-sections 
of the 1-phonon states. It is difficult to assess the effect of our > 3— phonon truncation 
on the calculation of the 2-phonon cross-sections. This will be investigated in a future 
publication. 

5.3.1 The Born approximation 

A striking disagreement between the predictions of the Born approximation and the exact 
calculation is displayed in Figure 9. At bombarding energy below about 2 A GeV, the Born 
approximation and exact calculation yield nearly the same cross-section for population of 
the 1-phonon M = state. However, at higher bombarding energy, the Born approx- 
imation cross-section becomes vanishingly small, whereas the exact calculation predicts 
appreciable cross-section. It must be recalled that the Born approximation involves only 
the on- shell Fourier component of the interaction, Vp a {u = h - ). It was shown in Ref- 
erence that AM = on-shell interaction matrix elements vanish at high bombarding 
energy in proportion to I/7 2 , as a result of cancellation between the contributions of the 
scalar and vector potentials. This accounts for the strong decrease at high bombarding 
energy of the Born AM = cross-section. However, off-shell interaction matrix elements, 
Vj3 a (uj 7^ E/3 ~ h Ea ), do not vanish at high 7; rather they diverge in proportion to log 7. In the 
exact calculation, off-shell interaction matrix elements contribute to the population of the 
one-phonon M = state, as do multistep processes. Thus the strong Born approximation 
selection rule against population of the 1-phonon M = state is not exhibited by the 
exact calculation. 

On the other hand, Figure 10 shows that the Born approximation and exact calculations 
are in fairly good agreement for the population of the 1-phonon \M\ — 1 state over the 
entire energy range. This suggests that most of the population of this state occurs in a 
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single-step process. The fact that the exact calculation predicts a slightly smaller cross- 
section can be interpreted as a result of the loss of flux from the \M\ = 1 state associated 
with the coupling to other available states. 

Figure 11 shows the cross-section for the total population of the 1- phonon level, 
including both the M = and the \M\ = 1 states. It is seen that as a result of the 
Born approximation underprediction for the M = state, and over-prediction for the 
\M\ — 1 state, there is fairly good agreement for the total 1-phonon cross-section over 
the entire energy range. Thus measurements of the total 1-phonon cross-section cannot 
distinguish between the predictions of the Born approximation and exact calculations, 
and our exact calculation has little to contribute to the many analyses of total 1-phonon 
cross-sections performed with the Born approximation. However, if the separate M = 
and \M\ = 1 cross-sections could be measured, then a clear choice could be made between 
the predictions of the Born approximation and the exact calculation. 

Figures 12 through 16 show our results for the 2-phonon states. In all cases, the 
cross-section predicted by the exact calculation far exceeds the predictions of the Born 
approximation. The natural explanation is that these states are predominantly excited by 
multi-step processes. 

5.3.2 Comparison with the harmonic vibrator model 

Some analyses of multiphonon excitations have used the harmonic vibrator model, in which 
it is assumed that the total effect of the interaction on the target protons can be imitated 
by an operator that is linear in the components of the collective variable R pn . This has 
the consequence that multi-phonon states are populated by a series of on-shell transitions, 
and thus the fact that the AM = on-shell transition matrix element becomes small at 
high bombarding energy implies that the AM = matrix elements will play little role 
in multiphonon excitation. It was shown in Reference [2] that the ratios of the cross- 
sections for the 2-phonon states (020), (021), (022), (100) are predicted by this model to 
be 1/0/3/2. The absence of population of the state (021)is due to the fact that this state 
consists of phonons with m = and m = 1, and the m = phonon is inaccessible. The 
predictions of the exact calculation, as shown in Figures 12 through 16, are quite different. 
At a bombarding energy of 10A GeV, the ratios are approximately 4.1/. 3/. 1/1.6. The 
important role of the AM = matrix elements in these numbers is evident from the 
bombarding-energy dependence shown in Figures 12 through 16. For example, of these 
four states, only (022) shows a cross-section that decreases with bombarding energy, and 
this is the state which could be expected to be least affected by AM = matrix elements. 
Thus, our microscopic treatment of the interaction, which follows the effect of the external 
pulse on each proton, gives quite different results from assumption that the interaction is 
linear in the collective variables. However, it must be borne in mind that the error we 
make by truncating our phonon space at two phonons may introduce significant errors 
into our 2-phonon predictions. This requires further investigation. 

5.3.3 The long-wavelength approximation (LWLA) 

It is seen in the cross-section comparisons, Figures 9 through 14, that at high bombarding 
energy the LWLA prediction exceeds the result of exact calculation, except for the 1- 
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phonon \M\ = 1 state and the 2- phonon \M\ = 2 state. The common characteristic of 
these two states is that they involve only oscillation in the transverse (x, y) directions, 
whereas all the other final states involve oscillation in the z direction. We have already 
noted in connection with Figures 4 and 6 that the excitation of the z oscillation takes place 
rapidly in the vicinity of t = 0. The LWLA matrix elements are more strongly peaked in 
the vicinity of t = (see Figures lb and 2b) due to absence of LWLA dynamic damping. 
Thus it is not surprising that the LWLA gives enhanced predictions for states involving 
z excitation. On the other hand, the excitation of the x,y oscillation is a more gradual 
process, covering a larger t-range, and so is more sensitive to the time region in which the 
exact matrix elements are greater than the LWLA matrix elements (Figures lb, 2b). 

6 Conclusions 

Although we have used a somewhat schematic representation of the GDR nuclear states, we 
can draw some conclusions about the extent to which the LWLA and Born approximation 
reproduce the results of our exact calculations. These conclusions may also apply to 
calculations that use more realistic nuclear states. 

It is apparent from Figures 9 through 16 that at bombarding energies below about 
3A GeV, the LWLA and the exact calculations predict essentially the same cross-sections, 
state by state. This observation supports the validity of several RCE investigations that 
have used the LWLA jHlEl- It is shown in Figure 12 that the agreement between the LWLA 
and exact calculations of the total 1-phonon cross-section extends up to about 5A GeV. 
Above 5A GeV, however, there are significant deviations between the LWLA and exact 
cross-sections. Thus at higher bombarding enrgyies, the LWLA ceases to be a reliable 
approximation, especially if it is important to know the cross-sections for population of 
individual M-states. 

Not surprisingly, the Born approximation is unable to predict the cross- section for 
population of 2-phonon states. Its performance with respect to the 1-phonon states is 
analagous to that of the LWLA. At bombarding energies below 2A GeV, it does well with 
both M = and \M\ = 1 states. At higher energies, it underpredicts the M = cross- 
section (Figure 9), and over-predicts the \M\ = 1 cross-section (Figure 10). The net result 
is good agreement with the exact prediction for the total 1-phonon cross-section up to 
about 8 A GeV (Figure 11). 
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Appendices 

A: Brink's representation of the GDR phonon states 

We first construct the states of the GDR which are automatically contained in the inde- 
pendent particle harmonic oscillator shell model. For simplicity, we restrict our attention 
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to N = Z closed-shell nuclei. We start by locating the Z protons and the Z neutrons 
relative to the fixed harmonic oscillator origin by p; , • • • , pz, , • • • , nz- To construct 
GDR states, we introduce new position variables 

R Pi H h Pz + nj j h n z 

2Z 

Rp „ = ( £ i ±_±P£) _ (=i ±^± !!£) (A1) 

and relative variables 7Tj , • • • , 7r^_i , f f , • • • , f z-j • The 7r j , • • • , 7r z _ J specify the lo- 
cation of the Z protons relative to their mass center. 




Figure Al. R pn locates the position of the proton mass-center relative to the neutron 
mass center, whereas R locates the total mass center. The j th proton is located relative 
to the proton mass center by fj. 



It is clear from Fig.(Al) that the jth proton is located relative to the origin by 

Pj = R + iR pn + f^TTj, •••,*■*-*) (A.2) 

We will not have to specify the tzx , ■ • • , ttz-i , Vi , • • • , vz-i or the fj any further; we 
will only need the R and the R pn dependence of pj. 

If the independent-particle harmonic oscillator Hamiltonian is expressed in terms of 
these new variables, it separates as follows: 

+ if' (7Tj, • • • ,7Tz-j, ^1, ' • ' ,V Z -l) • (A. 3) 
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Here m is the nucleoli mass and u sm is the single- nucleon oscillator frequency. Cor- 
responding to this separability, we can write our harmonic oscillator eigenstates in the 
form 



tt Q = *S' Ll (R) <f 2 ( R P n) x X^u-'-^z-u vu-'-^z-i) ■ (A.4) 



The spin degrees of freedom are incorporated in x, i n such a way as to produce a state which 
is antisymmetric with respect to the exchange of any two protons or any two neutrons. 
Note that changing N\, L\, M\ or N2, L2, M2 has no effect on the antisymmetry, since 
the variables R and R pn are symmetric with respect to the exchange of any two protons 
or any two neutrons. 

It is clear that the state of lowest energy, consistent with the Pauli principle, will have 
zero quanta in the variables R and R pn : 



Aground-state = *0 (R) *0 (Rpn) X X o{^l , ' ' ' , ™Z-1 , Vl,---,V Z -l) . (A.5) 



v^ L = <°(R) ^ L (R pn ) x X o{*i,---,*z-t, vi,---,"z-i) • (A.6) 



These states all have no quanta of center-of-mass motion. The relative motions of the 
protons to each other, and of the neutrons to each other, are exactly the same as they 
are in the shell-model ground state. The only feature distinguishing the states (A6) is 
the motion of the proton mass-center relative to the neutron mass-center. Brink [2H] 
identified these states as belonging to the GDR, because this collective oscillation of the 
protons relative to the neutrons was the original interpretation given by Goldhaber and 
Teller [57] to the phenomenon now known as the GDR. 

It can be seen from (A3) that the energy associated with the R pn quanta is huj sm , 
the same as the energy associated with the individual oscillator quanta. It was suggested 
by Wilkinson j2H], and proved by Brown and Bolsterli |31j . that the residual interaction 
between the nucleons will increase the GDR phonon energy. 

The harmonic oscillator radial functions associated with the individual nucleons were 
specified in (4.6). The same radial functions are used for the GDR states ^^ L (R pn ), 
except that m is replaced by Zm/2, which is the reduced mass for the relative motion of 
the proton and neutron spheres. 

Electromagnetic properties of states are normally discussed in the long-wavelength 
approximation, in which the electric dipole operator has the form 



Equation (A. 2) can be used to express this operator in terms of R, R pn , 7Tj. The only part 
of this operator that has matrix elements between states of the GDR band (A.6) is 



which can only connect states which differ by one GDR phonon. Thus the El transition 
between the ground state and the 1-phonon Brink level exhausts the El sum rule. 



Now let us consider a series of states of the form 





16 



B: Contribution of the scalar potential 

It is shown in j2Hj that if the nuclear states <f> a , cj)p are defined with "time- reversal" phases, 

(0kr = (-iy°- M « , (b.i) 

then the Vp a (uj) will be purely real. We will use this phase convention. 

The charge density pf3 a { Y ') is defined in terms of the nuclear eigenstates (ftp, 4> a by 

z 

p f3a (r')=e I d 3 Pl ...d 3 n z J2 S{r' - Pj )(f>*^ a (B.2) 

The sum extends over the protons only. By using (B2) and (4.3) in (4.1), we see that the 
scalar potential contribution to Vp a (uS) is 



|dV ^ et (r>) P/3a (r') = E<?A,»i^(^) fd 3 r' PPa {v') 3 ^r')Y^f') . (B.3) 



The integral required here also occurs in the Fourier transform of pp a {r') 



ppaQt) = J e lkr ' p Pa (r')d 3 r' = 4tt £ i A Yf(l) J d 3 r> p Pa (r')j x (kr')Y*(f') (B.4) 

A/x 

if we set k — \oj\/c. Thus we turn our attention to the evaluation of pp a (k): 

P/to (k) = J e lkr ' p Pa {v')d\' = J d\'^ r ' <^\eY, 8{v' -v j )\<f> a > 

i=l 

= <<Pp\ej2 e ik ^|0 a > . (B.5) 

We now choose states in the GDR band (A6) for the 4>p and <fi a , and use (A2) to express 
Pj in terms of the new variables: 

P^(k) = e j d 3 R\^°\ 2 X / d^n^^'^CRpO^lr^CRpn) 

x J d 3 7z u • • • , d 3 u z _ 1 \xo(^i, • • • ,nz-i, "1, • • • ,vz-i)\ 2 



z 

£ 

3=1 



Q y d 3 R pn ^'^(R pn ) e^^Z^iKn) (B.6) 



with 



Q = ej d 3 R\y° Q >°\ 2 e lkR x | dV,---,^-:*!^*, f,, 



1 2 



x V" e * f j(7ry, • • • j^-i, l/j, • • • , Vz-l) 

3=1 
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Note that the factor Q is independent of the GDR quantum numbers. This is a consequence 
of our assumption that the proton and neutron spheres remain unexcited during the R pn 
oscillation. Thus we can calculate Q by considering the ground-state proton density 



/1 - f°° 

d 3 rV k. r £ 2{2M)u 2 {rl) = 2 (2€+l) / u 2 n ,j(r')j (kry2dr' 

(B.8) 

The n',£ sum extends over the single-particle levels occupied in the ground state (for 
example: 00, 01, 02, 10, for 40 Ca). Using (B6) to evaluate poo(k) yields 



Poo(k) = e sz„ x Q , 



so that 



k 2 /-00 

Q = e^ xeJ2 2(2^ + 1) / u 2 n /r')j (kry 2 dr' . 
To evaluate the R pn integral in (B6), we write 

<^(R pn )*£^(R pn ) = (-i)^-^ ^(R pn )K a : La (^) 



(B.9) 
(B.10) 



= ( _ 1) l,-m, ^ _ Mp M a \LM a -Mp)x [* Jv ^ L ' 3 (R pn )* JV - i «(R pn ) 

L 



L 

M a —M{3 ' 

(B.H) 



We then use the Talmi-Moshinsky coefficients, defined by 

L 



^ ni£l (n)^ 2 (r 2 



"2^2 / 



M 



E « 

n,\,N,A 

,ri - r 2 



^2 



ni£i,n 2 £ 2 |nA, iVA) L 
) 



V2 



M 



;b.i2) 



with ri = r 2 = R pn , to write 

^MRpn)^' L "(R P n) 



L 



= £ ^ q +l,-a-a (Ar /3 L /3 ,iV a L a |nA,iVA) L [ ^" A (0) ^(v^R, 

n,X,N,A 



pnj 



L 

M a -Ma 



^ 2 " L (^) 3/4 E (A^,A^>0,^ (B.13) 



ra,JV 



The extra factor ^i+^-a-a j s ^ ue ^ our uge Q f time-reversal phases. The n, A/" sum in 
(B13) is restricted by the quanta-preserving condition 2(n+N)+L = 2(N /3 +N a ) + Lp+L a , 
which in turn requires that Lp + L a — L be even. 

To Fourier transform ^j^_ M (-\/2Rp n ) we use the general harmonic oscillator result 



1 /■« 

^ y_ c 



(2tt) 3 / 2 



e lk ' r ^(r)rf 3 r 



•2n+^ „/,ri£ 



C(k) , 
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where ip 1 ^ f(k) has a size parameter v which is the reciprocal of the size parameter of ip™ (r). 
Then 

jT e^- R -^_ M/3(y 2R pn)rf 3 Rpn = ^ ^^^^ 

= (-1)^ tt 3 / 2 y^_ M/3 (k)fe(^). (B.i4) 



By combining (B6), (B8), (B13) and (B14), we get 

n',£ 

x ^ i 2N + L (LpL a — MpM a \L M a — Mp) 



pp a (k) = (-l)^^-M^) 3/4 e^ E 2 (2^ + l) r jo(^'K v (r'y 2 rfr' 



x(N $ L p ,N a L a \na,NL) L x sj ^J^ u NL {^j=)Y^_ M ^) (B.15) 
Finally, comparison with (B3) and (B4) shows that 

r j3 I ret// \ //x {-l) M H La ~ L P ,Zv7T. 3/4 (k/c) 2 ( 

y dV ^ c ct (r',^)p /3a (r / ) = ^ (—) ' e 8Z " ^(^^7) 

/•oo \u}\ 

x £2(2^ + 1) / j ( U r') u 2 n/ (r'y 2 dr' ]T C l , M/3 -m q — MpM a \L M a — Mp) 



- 

X 

n,7V 



l(9n 4- lTn H 

£(-1)" (A^, A a L>0,7VL) L ^ J,, J - ^l(^), (B.16) 



which gives us an explicit expression for the contribution of the scalar potential to Vp a (cu). 
This expression is summarized in (4.5) 

C: Contribution of the vector potential 

The analysis is similar to the one presented in Appendix B for the scalar potential. We 
need 

-\ /rfVj^(r')-Af(r',uO = ~ j d 3 r'\jp a (v%^(v',u;) 

= -JrE Ca,, / d 3 r h(kr')Y^v%p a (v% .(C.l) 

To obtain these integrals we study the Fourier transform of [j/3 a (r')] z : 



[j/Ja(k)], = / e^ r ' [jp a (r% d 3 r' = Att^ *%* A (k) f <Pr j x (kr>)\j fla (i')], 

= / rfV."'<i 3 «z E e ,k ' PJ x U' -^-K - <t> a J-<l>"e) (C2) 
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We now take cf>p, <p a to be states in the GDR band (A6), and transform to the variables 
R, R pn , izi , ■ • • , 7Zz-i ? , • • • , vz-i • To do this we use (A2) to obtain 







1 d 



dp 



3,z 



2ZdR 7 



+ 







ZdR, 



+ 



E 



k=l^j-l;n=x,y,z 



( 9 







Since the it/- involve only the relative positions of the Pj, the quantities 
independent of R and R pn 

Equation (C2) now takes the form 



-n k ^ are 



[j/?a( k ) 



eh 



2(Zmi) 



Ox/ 



X 



5 



•9 j*N a L f 



(Rpn) ) (C3) 



L pn,z <V - I, 'pn,2 

where Q is the quantity defined in (B7) and evaluated in (BIO). Note that (C3) vanishes 
for diagonal matrix elements (Np = N a , Lp = L a , Mp = M a ). The derivatives in (C3) 
can be expressed as linear combinations of harmonic oscillator wave functions using the 
relation 



d_ 

dZ 



X 



(e+l-m)(e + l + m) 
\ (2£+l)(2£ + 3) 



2n^"" 1 ' m (R) + v / 2n + 2£ + 3^ £+1 (R) 



This converts (C3) into a linear combination of terms such as (B6), which can be evaluated 
exactly as was done in Appendix B. 



D: Symmetries 

Since the projectile moves in the y — z plane, its scalar and vector potentials will be 
invariant under reflection across the y — z plane. Therefore, only reflection-symmetric 
target states will be excited from the reflection-symmetric ground state. These are 

\b NL or 1 (\b NL 4- \b NL \ 

> or "^2 |M| ~\ M \>' 

Thus if M a ^ and Mp ^ 0, we must calculate 

< + ^/)|VH|i=(^ + > 



> 
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If M/3 ^ and M a = 0, we must calculate 

It can be verified that Equation B3 of Reference 

<^\v(u)\*»i L «> = <Kl La \v(-u)\*Zf">, 

is satisfied by the matrix elements discussed in Appendices B and C. 
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Figure Captions 

Figure 1. (a) Vooo-on(^) for a bombarding energy of 10A GeV and an impact pa- 
rameter of 12 fm. The peak at oo = is narrower for the exact curve than for the long- 
wavelength approximation because the exact curve shows the effect of dynamic damping, 
(b) Vooo-oii(^) corresponding to the Vooo-onO^) of Figure la. The narrower u dependence 
of the exact curve results in a broader t dependence. 

Figure 2. (a) The same as Figure la, but for an impact parameter of 50 fm. (b) The 
same as Figure lb, but for an impact parameter of 50 fm. 

Figure 3. Excitation probability of the state ^[^/5 1 + ^-\] as a function of t, for a bom- 
barding energy of 10A GeV and an impact parameter of 200 fm. The Born approximation 
excitation probability is shown for comparison. 

Figure 4. The same as Figure 3, except that the state is ^q 1 . The t — > oo prediction 
is 0.6 x 10 -5 , compared to a Born approximation prediction of about 10~ 8 . 

Figure 5. The same as Figure 3, but for an impact parameter of 12 fm. 
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Figure 6. The same as Figure 4, but for an impact parameter of 12 fm. 

Figure 7. The asymptotic excitation probability of the state ^[^? 1 + x I /01 i] as a function 
of impact parameter, for a bombarding energy of 10A GeV. 

Figure 8. The same as Figure 7, except that the state is -j= [^Sp + \I/° 2 2 ] ■ Because of the 
importance of multistep processes in the excitation of this state, the Born approximation 
underestimates the excitation probability even at large impact parameter. 

Figure 9. Cross-section for excitation of the state ^q 1 as a function of the bombarding 
energy. 

Figure 10. Cross-section for excitation of the state -j= [^f^ 1 + ^°\] as a function of the 
bombarding energy. 

Figure 11. Total cross-section for excitation of the one-phonon states (the sums of the 
curves in Figures 9 and 10). 

Figure 12. Cross-section for excitation of the state as a function of the bombarding 
energy. The maximum Born cross-section is 0.038 mb. 

Figure 13. Cross-section for excitation of the state ^75 [^i 2 + ^ 02 i] as a function of the 
bombarding energy. 

Figure 14. Cross-section for excitation of the state ^75 + ^- 2 2] as a function of the 
bombarding energy. 

Figure 15. Cross-section for excitation of the state \1/q as a function of the bombarding 
energy. The maximum Born cross-section is 0.0006 mb. 

Figure 16. Total cross-section for excitation of the two-phonon states (the sums of the 
curves in Figures 12, 13, 14, and 15). The maximum Born cross-section is 0.65 mb. 
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